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1 3.  ABSTRACT  ( Maximum  200  words) 

Hyperspectral  imaging  cameras  require  smooth  motion  in  order  to  acquire  uncorrupted  images.  Typically,  a  camera  is  mounted  on  a 
satellite  that  then  rotates  to  allow  the  camera  to  scan  the  earth’s  surface.  The  rotation  of  the  satellite  causes  propellant  sloshing  in  its  fuel 
tank,  which  can  disrupt  the  motion  of  the  satellite  and  therefore  degrade  the  quality  of  the  hyperspectral  images. 

The  severity  of  the  liquid  sloshing  has  been  investigated  with  a  numerical  fluid  dynamics  model  called  Flow-3 D  (Flow  Science, 
1997).  This  model  is  well  suited  for  such  a  problem  and  has  been  previously  verified  by  comparing  numerical  predictions  with  experi¬ 
mental  slosh  data  in  reduced  gravity  environments  (Sicilian  and  Tegart,  1989,  Fisher  et  al.,  1992).  The  experimental  data  were  collected 
from  drop  tower  tests,  which  produce  a  controlled  acceleration  environment. 
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NUMERICAL  SIMULATION  OF  LIQUID  SLOSH  IN  MICROGRAVITY 


INTRODUCTION 

Hyperspectral  imaging  cameras  require  smooth  motion  in  order  to  acquire  uncorrupted 
images.  Typically,  a  camera  is  mounted  on  a  satellite  that  then  rotates  to  allow  the 
camera  to  scan  the  earth’s  surface.  The  rotation  of  the  satellite  causes  propellant  sloshing 
in  its  fuel  tank,  which  can  disrupt  the  motion  of  the  satellite  and  therefore  degrade  the 
quality  of  the  hyperspectral  images. 

The  severity  of  the  liquid  sloshing  has  been  investigated  with  a  numerical  fluid  dynamics 
model  called  Flow-3D  (Flow  Science,  1997).  This  model  is  well  suited  for  such  a 
problem  and  has  been  previously  verified  by  comparing  numerical  predictions  with 
experimental  slosh  data  in  reduced  gravity  environments  (Sicilian  and  Tegart,  1989, 
Fisher  et  al.,  1992).  The  experimental  data  were  collected  from  drop  tower  tests,  which 
produce  a  controlled  acceleration  environment. 

PROBLEM  DESCRIPTION 

In  the  numerical  simulations,  the  satellite’s  fuel  tank  was  the  domain  of  interest.  The 
simulations  began  with  an  initial  liquid  configuration  and  the  tank  was  then  accelerated  in 
a  prescribed  manner  corresponding  to  a  typical  satellite  rotating  motion.  The  sloshing 
fluid  exerts  forces  on  the  tank,  which  can  be  computed  by  the  numerical  model  and 
assessed  for  their  impact  on  the  satellite’s  motion. 

Two  different  masses  of  liquid  have  been  considered  in  order  to  investigate  the  influence 
of  fuel  volume  on  the  resulting  forces  acting  on  the  tank.  The  first  case  is  a  high  mass 
case  with  59.5  kg  of  propellant,  which  corresponds  to  a  liquid  volume  of  0.0593  m3.  The 
remaining  volume  of  the  tank  is  filled  with  gas  at  a  pressure  equal  to  1.29  MPa.  Fig.  1 
shows  three  slices  of  the  geometry  of  the  tank  along  with  the  initial  fluid  configuration. 
Note  that  z  denotes  the  vertical  direction,  x  andy  denote  the  two  lateral  directions,  and  the 
origin  is  at  the  satellite’s  center  of  mass.  The  red  area  represents  the  propellant,  while  the 
blue  area  is  gas.  The  tank  is  spherical  in  shape  with  a  diameter  equal  to  0.5588  m  and 
there  is  a  cylindrical  sump  at  the  bottom  with  a  diameter  of  0.2032  m  and  height  of  0.127 
m.  There  are  also  two  sets  of  vanes,  which  are  perpendicular  to  each  other,  running 
vertically  around  the  circumference  of  the  tank.  The  vanes  are  approximately  0.01905  m 
high.  The  second  case  is  a  low  mass  case  with  30  kg  of  fuel,  which  corresponds  to  a 
volume  of  0.0299  m3  and  the  gas  pressure  is  783  kPa.  This  case  is  illustrated  in  Fig.  2. 

In  a  low  gravity  environment,  surface  tension  is  the  dominant  restoring  force  acting  on 
the  liquid-gas  interface.  This  is  expressed  in  terms  of  the  Bond  number,  which  is  the 
ratio  of  the  gravity  to  surface  tension  forces  and  is  written  as 


Bo  = 


pgR 2 


(1) 
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where  p  is  the  liquid  density,  g  is  the  gravitational  acceleration,  R  is  the  radius  of  the 
tank:,  and  <y  is  the  surface  tension  coefficient.  The  liquid  is  modeled  as  anhydrous 
hydrazine  and  has  a  density,  p =1004.0  kg/m3  and  dynamic  viscosity,  p  =  0.0009  kg/m  s 
at  25  degrees  C,  while  <7  =0.006667  kg/s2.  In  addition,  the  wall  contact  angle  is  40 
degrees.  In  the  present  simulations,  g=0  was  specified  and  therefore  Bo= 0,  which 
indicates  the  dominant  role  of  surface  tension.  Therefore,  surface  tension  was  included  in 
the  numerical  simulations. 

The  Reynolds  number  is  another  important  parameter  in  fluid  dynamics.  It  is  the  ratio  of 
inertia  to  viscous  forces  and  can  be  written  as 


(2) 


where  u  is  a  typical  liquid  velocity,  L  is  a  length  scale  of  the  flow,  and  v=p/p  is  the 
kinematic  viscosity  of  the  liquid.  The  Reynolds  number  also  indicates  the  importance  of 
turbulence.  A  flow  with  a  larger  Re  is  more  turbulent  than  a  similar  flow  with  a  smaller 
Re.  The  numerical  simulations  revealed  that  typical  fluid  velocities  were  on  the  order  of 
0.005  m/s,  while  a  typical  length  scale,  taken  as  the  depth  of  flow  in  the  tank,  was 
approximately  0.3  m.  This  yields  an  Re  of  about  1500,  which  indicates  that  the  flow  is 
laminar.  Therefore,  turbulence  was  not  considered  in  the  simulations. 

THE  GOVERNING  EQUATIONS 

The  liquid  motion  is  modeled  with  the  conservation  of  momentum  equations  in  three 
orthogonal  directions.  In  addition,  the  liquid  is  assumed  to  be  incompressible,  and 
therefore  its  density  remains  constant.  Also,  as  previously  discussed,  the  flow  is  assumed 
to  be  laminar.  Under  these  assumptions,  the  momentum  conservation  laws  reduce  to  the 
so-called  Navier-Stokes  equations,  which  are  written  as 


du 

du 

du 

du 

+  u — 

+  v — 

+  w — 

dt 

dx 

dy 

dz 

-1  dp 
p  dx 


+  v 


f  d2u  d2u  d2u 
dx 2  dy 2  dz2 


+  ar 


(3) 


dv  dv  dv 

—  +  u — +v - b  w- 

dt  dx  dy  dz 


dv  -l  dp  ( d2v  d2v  92v^ 
—  = —  +  V\  — T  +  ^  +  T-T 


p  dy 


dx  dy  dz 


+  a„ 


dw  dw  dw  dw  _  -1  dp  ( d2w  d2w 

&r+lF+aF' 


-  +  U - +  V - +  W - = - - — bV\ 

dt  dx  dy  dz  p  dz 


(4) 

(5) 


where  u,  v,  and  w  are  the  velocities  in  the  x,  y,  and  z  directions,  respectively,  while  p 
represents  the  pressure.  Note  that  the  equations  are  written  in  terms  of  a  coordinate 
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system  that  moves  with  the  tank.  The  terms  a„  ay,  and  a2  denote  the  components  of  the 
acceleration  of  the  tank  reference  frame  with  respect  to  an  inertial  reference  frame.  If  the 
moving  reference  frame  undergoes  angular  accelerations  about  its  origin  only,  the 
acceleration  terms  become 


ax  =  -| <oy  z-o)zy  +  2(& y w  -  o)zv)  +  a)y  (axy  -  a>yx)  -  mz  (azx  -  axz)j  (6) 

ay  =  -|  G)z  x  —  co x  z  +  2(g)zu  -  c oxw )  +  az(coyz  -  (azy)  -  cox(coxy  -  coyx)  |  (7) 

az  =-\o)xy-G)yx  +  2(coxv  -coyu)  +  cox  (cozx  -coxz)-  coy  (o)yz  -  G)zy)\  (8) 


where  a  denotes  the  angular  velocity  and  the  dot  denotes  the  angular  acceleration.  In  the 
present  simulations,  only  the  angular  velocity  and  acceleration  about  the  y  axis  are 
nonzero.  The  specific  form  of  Oy  is  shown  in  Fig.  3,  while  the  angular  acceleration  is 
calculated  by  numerically  differentiating  the  angular  velocity  with  respect  to  time. 

The  conservation  of  liquid  mass,  under  the  incompressible  flow  assumption,  yields  the 
additional  constraint, 


du  dv  dw  .  ('9') 

dx  dy  dz 


Eqs.  (3)  through  (5)  and  Eq.  (9)  form  a  coupled  set  of  nonlinear  partial  differential 
equations  to  be  solved  for  the  unknowns  u,  v,  w,  and  p  as  functions  of  space  and  time. 
The  resulting  shear  stresses  and  pressure  can  then  be  integrated  over  the  wall  of  the  tank 
in  order  to  estimate  the  resulting  net  forces  and  moments  induced  by  the  sloshing  fuel. 


The  movement  of  the  liquid-gas  interface  is  tracked  using  an  Eulerian  approach  by 
measuring  fluid  movement  in  and  out  of  the  stationary  computational  cells.  This  method 
was  first  developed  by  Hirt  and  Nichols  (1979)  and  is  known  as  the  Volume  of  Fluid 
(VOF)  method.  The  fluid  movement  is  tracked  by  solving  the  following  transport 
equation, 


dF 

dF 

dF 

dF  . 

+  K - 

+v — 

+  w —  =  0 

dt 

dx 

dy 

dz 

(10) 


where  F  is  the  ratio  of  the  liquid  volume  in  a  given  cell  to  the  total  volume  of  that  cell. 
This  approach  can  model  discontinuous  interfaces  as  well  as  liquid  detachment,  which 
can  be  extremely  important  in  the  simulation  of  sloshing  fluids.  Note  that  this  approach 
neglects  the  motion  of  the  gas  and  therefore  implicitly  assumes  that  the  gas  motion  has  a 
negligible  effect  on  the  movement  of  the  fuel.  Also,  the  gas  pressure  is  assumed  to 
remain  constant  during  the  simulation. 
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NUMERICAL  SOLUTION 


The  governing  equations  cannot  be  analytically  solved  for  a  complex  problem  involving 
sloshing  fluids.  Therefore,  the  equations  must  be  solved  numerically.  As  previously 
mentioned,  the  numerical  model  employed  in  this  study  is  commercially  available  and  is 
called  Flow-3D.  It  is  a  finite-difference  model  that  is  based  on  the  Marker- And-Cell 
method  developed  by  Harlow  and  Welch  (1965).  All  temporal  and  spatial  derivatives  in 
Eqs.  (3)  through  (5),  Eq.  (9),  and  Eq.  (10)  are  replaced  with  finite  difference 
approximations  in  order  to  obtain  a  discrete  solution.  The  model  uses  a  stationary, 
Cartesian  computational  grid  that  may  be  selectively  refined  in  desired  areas.  The  fuel 
tank  geometry  is  incorporated  into  the  computational  domain  by  linearly  blocking  out 
portions  of  cells  that  fall  outside  of  the  walls  of  the  tank.  Corresponding  values  of  cell 
volumes  and  faces  that  are  open  to  flow  are  then  computed  and  used  in  the  discrete 
equations  to  compute  the  spatial  gradients  of  velocity,  pressure,  and  F  (Hirt  and  Sicilian, 
1985).  The  grid  is  staggered,  with p  and  F  defined  at  cell  centers,  while  the  velocities  are 
defined  at  cell  faces,  as  shown  in  Fig.  4.  This  prevents  “checkerboard”  solutions  that  can 
result  in  incompressible  flow  computations  when  the  velocity  and  pressure  become 
unlinked  at  adjacent  computational  nodes  (Patankar,  1980). 

The  model  solves  the  momentum  and  incompressibility  equations  in  an  iterative  manner 
using  a  pressure  correction  scheme.  A  computationally  efficient  scheme  may  be  obtained 
by  differentiating  Eqs.  (3)  through  (5)  in  space  and  substituting  the  results  into  Eq.  (9). 
This  yields  a  Poisson  equation  for  p  at  the  new  time  level  n+1.  This  equation  is  solved 
by  Successive  Overrelaxation  (SOR)  and  the  new  pressure  is  then  used  to  compute  new 
velocity  estimates  by  solving  the  discrete  analogs  of  Eqs.  (3)  through  (5).  The  new 
velocity  estimates  are  then  used  to  compute  a  new  pressure  field  via  the  Poisson  equation 
and  this  process  is  continued  until  the  change  in  p  becomes  smaller  than  some  tolerance. 
At  this  point,  the  new  flow  field  is  known  and  is  then  used  to  compute  the  new  free 
surface  configuration  by  solving  the  discrete  version  of  Eq.  (10).  Boundary  conditions 
are  then  applied  to  all  flow  parameters,  which  completes  the  solution  for  one  time  step. 

Boundary  conditions  are  applied  as  follows.  At  the  walls  of  the  tank,  the  velocity 
perpendicular  to  the  wall  is  set  to  zero.  Velocities  tangent  to  the  wall  are  set  according  to 
the  local  wall  shear  stress,  which  depends  on  the  local  Re.  At  the  liquid-gas  interface, 
velocities  outside  the  liquid  are  extrapolated  from  the  interior  such  that  the  resulting 
tangential  shear  stress  is  zero.  The  pressure  in  a  cell  containing  a  liquid-gas  interface  is 
extrapolated  from  the  interface.  Surface  tension  is  converted  to  an  equivalent  pressure 
that  is  specified  at  the  liquid-gas  interface.  When  a  solid  wall  is  in  contact  with  a  liquid- 
gas  interface,  the  wall  adhesion  force  is  also  estimated  using  the  specified  wall  contact 
angle.  The  net  force  is  then  converted  to  an  equivalent  interface  pressure. 

The  domain  was  divided  into  125,000  equal-sized,  cube-shaped  cells.  The  grid  spacing 
in  the  x  direction,  Ax,  was  0.01 1  m.  The  grid  spacings  in  the  other  directions,  Ay  and  Az, 
were  the  same.  The  time  increment,  At,  varied  during  the  computation  and  was 
automatically  chosen  by  the  model  to  maintain  the  stability  and  accuracy  of  the  solution. 
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Additional  simulations  were  performed  with  216,000  cells  in  order  to  check  the 
convergence  of  the  model.  Also,  more  accurate  discrete  approximations  to  the  spatial 
and  temporal  derivatives  were  employed  as  an  additional  test.  Neither  refinement  led  to 
significant  changes  in  the  numerical  solutions. 

MODEL  RESULTS 

Fig.  5  shows  plots  of  the  pressure  contours  and  velocity  vectors  for  the  high  fuel  mass 
case.  The  plots  are  for  the  x-z  plane  at  y= 0.  The  first  frame  shows  the  flow  field  at  *= 40 
s.  At  this  point,  the  tank  has  just  finished  moving  upward  and  has  just  begun  its  descent. 
It  is  seen  that  the  fuel  is  continuing  to  move  in  a  counterclockwise  direction  due  to  inertia 
and  much  of  the  fuel  has  moved  away  from  the  top  of  the  tank  down  toward  the  sump. 
The  next  frame  is  at  *=96  s,  at  which  time  the  tank  has  just  finished  accelerating 
downward,  but  continues  downward  at  a  constant  angular  velocity.  The  fuel  has  been 
dragged  by  the  tank  walls  into  a  clockwise  motion,  and  fuel  has  moved  back  toward  the 
right  side  of  the  tank.  At  *= 228  s,  the  downward  movement  has  been  completed  and  the 
tank  is  accelerating  upward  again.  Again,  the  fuel  motion  is  still  clockwise  due  to  inertia. 
At  *=280  s,  the  tank  has  completed  its  upward  motion  and  is  beginning  to  descend  again. 
A  strong  counterclockwise  flow  induced  by  the  tank  motion  is  clearly  evident  in  the 
figure.  At  *=320  s,  the  tank  has  completed  its  downward  motion  and  has  come  to  rest. 
However,  the  strong  counterclockwise  current  has  only  been  damped  by  the  downward 
movement  of  the  tank.  At  *=400s,  there  is  still  significant  fuel  motion.  Slowly,  viscous 
damping  would  dissipate  the  remaining  kinetic  energy  of  the  liquid  and  surface  tension 
would  restore  the  liquid-gas  interface  to  its  equilibrium  configuration. 

Fig.  6  shows  plots  of  the  pressure  contours  and  velocity  vectors  for  the  low  fuel  mass 
case.  These  plots  are  for  the  same  plane  as  the  previous  case.  The  fuel  behaves  in  much 
the  same  manner  as  in  the  previous  case.  However,  the  smaller  fuel  volume  allows  some 
areas  of  the  wall  of  the  tank  to  become  dry. 

Fig.  7  compares  the  forces  acting  on  the  fuel  tank  in  the  x,  y,  and  z  directions  for  the  two 
cases.  Fx  denotes  the  force  in  the  x  direction,  while  Fy  and  Fz  denote  the  forces  acting  in 
the  y  and  z  directions,  respectively.  As  expected,  Fx  and  Fy  are  quite  small  and  have  no 
discemable  harmonic  components  since  most  of  the  fluid  motion  was  in  the  z  direction. 
Fz  is  slightly  larger  in  magnitude  and  there  is  also  a  more  distinct  harmonic  nature  to  the 
force,  which  is  attributable  to  the  satellite’s  motion.  There  is  no  fundamental  frequency 
since  the  only  restoring  force  in  this  environment,  surface  tension,  is  quite  small.  Fig.  8 
shows  the  moments  about  the  origin  in  the  x,  y,  and  z  directions.  Mx  and  Mz  are  somewhat 
smaller  than  My,  which  is  again  expected  because  the  tank  was  rotated  about  the  y  axis. 


SUMMARY  AND  CONCLUSIONS 

A  computational  fluid  dynamics  model  has  been  used  to  simulate  the  sloshing  of 
propellant  in  a  satellite  fuel  tank  in  a  microgravity  environment.  The  resulting  forces 
acting  on  the  tank  were  computed  in  order  to  evaluate  their  impact  on  the  motion  of  the 
satellite.  The  sloshing  was  created  by  a  rotating  motion,  which  was  chosen  to  simulate 
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the  scanning  motion  of  the  satellite  needed  for  collecting  hyperspectral  images.  A 
commercial  code  called  Flow-3  D  was  utilized  to  perform  the  simulations.  The  code  has 
been  previously  verified  for  sloshing  liquid  problems  in  reduced  gravity  environments 
and  therefore  no  verification  was  performed  in  this  study.  However,  numerical 
experiments  with  additional  grid  refinement  and  more  accurate  discrete  approximations 
of  the  derivatives  were  performed  to  verify  that  the  solutions  presented  here  were 
converged. 

In  general,  the  resulting  forces  and  moments  acting  on  the  tank  were  quite  small  and 
should  have  no  adverse  impact  on  the  motion  of  the  satellite.  The  harmonic  nature  of  Fz 
and  My  was  due  to  the  satellite  motion,  and  not  the  restoring  force  of  surface  tension.  In 
addition,  the  amount  of  fuel  in  the  tank  had  little  impact  on  the  magnitudes  of  the  forces 
and  moments.  These  results  are  expected  due  to  the  rather  gentle  motion  of  the  fuel  tank. 
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APPENDIX 

The  previous  simulations  were  performed  in  an  environment  with  no  gravity  acting  on 
the  tank  or  liquid.  Therefore,  surface  tension  was  the  dominant  force  that  acted  to  restore 
the  system  to  equilibrium  conditions.  Since  this  force  is  quite  weak,  no  dominant 
harmonic  character  was  observed  in  the  slosh  forces  acting  on  the  tank.  In  order  to  better 
demonstrate  this,  three  additional  simulations  were  performed  for  the  same  tank 
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undergoing  a  constant  acceleration  in  the  vertical  (z)  direction.  Three  different 
accelerations  equal  to  0.00068g,  g,  and  3 g  were  utilized.  The  tank  was  then  impulsively 
perturbed  in  the  x  direction  and  the  resulting  slosh  forces  were  computed.  The  power 
spectral  densities  (PSD)  of  the  x  component  of  the  slosh  force  are  plotted  in  Fig  9.  As  a 
function  of  frequency.  This  figure  shows  that  the  low  acceleration  case  has  a  peak  in  the 
PSD  at  a  lower  frequency  than  the  g  acceleration  case.  Likewise,  the  g  case  has  its  peak 
PSD  value  at  a  lower  frequency  than  the  3 g  case.  This  is  because  the  larger  acceleration 
creates  a  more  dominant  restoring  force,  which  acts  more  rapidly  in  bringing  the  system 
to  equilibrium.  This  in  turn  leads  to  more  rapid  fuel  slosh  and  therefore  a  higher 
frequency  in  the  peak  of  the  PSD. 
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Figure  1:  Initial  fuel  configuration  for  the  high  mass  case  (units  in  meters). 
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Figure  5:  Pressure  contours  and  velocity  vectors  for  the  high  fuel  mass  case 
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Figure  6:  Pressure  contours  (Pa)  and  velocity  vectors  (m/s)  for  the  low  mass  case. 
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Power  Spectral  Density 


Frequency  (Hz) 

Figure  9:  Power  spectral  density  of  the  x  slosh  force. 
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